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Definition 

Normalized-Difference Snow Index (NDSI) - normalized 
difference of two bands (one in the visible and one in the 
near-infrared or short-wave infrared parts of the spectrum) 
is used to map snow. Snow is highly reflective in the visi- 
ble part of the EM spectrum and highly absorptive in the 
near-infrared or short-wave infrared part of the spectrum, 
whereas the reflectance of most clouds remains high in 
those same parts of the spectrum, allowing good separa- 
tion of most clouds and snow. »! If 


Introduction 


The NDSI has a long history. The use of ratioing visible 
(VIS) and near-infrared (NIR) or short-wave infrared 
(SWIR) channels to separate snow and clouds was 
documented in the literature beginning in the mid-1970s 
by Valovcin (1976, 1978) and also by Kyle et al. (1978). 
A considerable amount of work on this subject was 
conducted at, and published by, the Air Force Geophysics 
Laboratory (AFGL) (e.g., see Bunting and d'Entremont, 
1982), The objective of the AFGL work was to discrimi- 
nate snow cover from cloud cover using an automated 
algorithm to improve global cloud analyses. Later, auto- 
mated methods that relied on the V1S/N1R ratio were 
refined substantially using satellite data, by Crane and 
Anderson (1984), Dozier (1989), and Rosenthal and 
Dozier (1996) for regional scales, and by Riggs et al. 



(1993), Hall et al. (1995, 2002), and Hall and Riggs 35 
(2007) for global snow-cover mapping. In this section, 36 
we provide a brief history of the use of the NDSI for map- 37 
ping snow cover. 38 

Band ratios used to discriminate snow and clouds 39 

Results of an investigation of snow reflectance character- 40 
istics using data from Skylab Earth Resources Experiment 4i 
Package (EREP) S 1 92 multispectra! scanner are presented 42 
by Barnes and Smallwood (1975). For the first time, satel- 43 
lite study of snow from the spectral range extending from 44 
the VIS to the 1R (0.41-12.5 pm) was possible, and this 45 
paved the way for automated snow-cover mapping. 46 
Shortly thereafter, Valovcin (1976) at AFGL introduced 47 
the idea of using the ratio of radiance values in the VN1R 48 
(0.68-0.76 pm) and NIR or SWIR (1 .55 1 .75 pm) to pro- 49 
vide a method to discriminate between snow' cover and 50 
clouds. Kyle et al. (1978) used the ratio of the 1.6- 51 
0.754 pm channels to distinguish snow and clouds using 52 
a cloud physics radiometer with 0.754 1.64 pm channels. 53 
They also used an 1R band to test for surface temperature 54 
further distinguished snow and clouds. 55 

Additional work done at AFGL by Bunting and 56 
d’Entremont (1982) employed a 1.6 pm sensor flown on 57 
the Defense Meteorological Satellite Program (DMSP) 58 
Special Sensor C (SSC) to separate snow and clouds. They 59 
also used 11% reflectance to define the lower bound of 60 
reflectance for snow cover. Crane and Anderson (1984) 6i 
reviewed the previous work, mainly conducted at AFGL, 62 
and employed the DMSP Operational Linescan System 63 
(OLS), which operated in the 0.4 1.0 pin and 8-13 pm 64 
range, along with SSC data (1.51-1.63 pm). They 65 
employed reflectances derived from the various sensors 66 
to map snow using a threshold technique. 67 

More-sophisticated use of band ratios as applied with 68 
Landsat Thematic Mapper TM data was developed by 69 
Dozier (1987, 1989). The nonnalized difference of TM 70 
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71 bands 2 (0.52-0.60 |tm) and 5 (1.55 1.75 |im) was intro- 

72 duced in Dozier (1989). Dozier and Marks (1987) discuss 

73 automated snow mapping and threshold tests for 

74 shadowed snow, cloud, vegetation, and soil in sunlit areas. 

75 With the anticipated launch of the Moderate Resolution 

76 Imaging Spectroradiometer (MODIS) at the end of the 

77 1990s, a global snow-mapping algorithm needed to be 

78 developed that would perform automatically and not be 

79 computationally intensive. Using the heritage algorithms 
so discussed above, Hall et al. (1995) coined the term 

81 normalized-difference snow index and outlined a snow- 

82 mapping algorithm that would be the basis of the MOD1S 

83 standard snow-mapping product. The prototype algo- 

84 ritbm, called Snowmap, used a normalized difference 

85 between MOD1S band 4 (5.45-5.65 pm) and 6 (1.628- 

86 1.652 pm), as was done in Bunting and d'Entremont 

87 (1982), Crane and Anderson (1984), and Dozier (1989) 

88 using TM bands 2 and 5. The prototype MODIS algorithm 

89 also employed several spectral tests. A planetary reflec- 

90 tance < 1 1 % was a threshold test in which values <11% 

91 were mapped as “not snow,” determined not to be snow. 

92 The prototype MODIS snow-mapping algorithm was 

93 improved with additional spectral tests. One key modifica- 

94 tion is that the NDSI threshold was changed in forested 

95 areas based on results of a canopy reflectance model 

96 (Klein et al., 1 998), using both the Normalized Difference 

97 Vegetation Index (NDVI) and NDSI in densely forested 

98 areas as determined from the NDVI test. A thermal mask 

99 was also included to remove erroneous “snow” in loca- 

100 tions where snow is considered to be impossible. Small 

101 specks of erroneous snow that show up on an image may 

102 be due to sand. If the band 31 (10.780-1 1 .280 pm) tem- 

103 perature is >283 K, then a pixel is considered “not snow;” 

104 This type of thermal test of surface temperature had previ- 

105 ously been used by Kyle et al. (1978) and Romanov and 

106 Gutman (2000). The standard MODIS cloud mask is also 

107 employed as an input to the snow algorithm. 

106 Following the 1999 launch of the MODIS on the Terra 

109 spacecraft, the snow algorithm was modified several 

110 times, but the NDSI has remained the basis of the algo- 

111 rithm. The current algorithm is Version 005 (see Riggs 

112 etal., 2006). 

113 Summary 

114 The term normalized-difference snow index (NDSI) was 

115 coined by Hall et al. (1995), but the NDSI technique 
lie already had nearly a 20-year heritage as similar methods 

117 using various visible and near-infrared bands had been 

1 18 used since the mid- 1 970s to map snow and separate snow 

119 from most clouds. Following the launch of the MODIS in 

120 1 999, the NDSI approach 10 mapping snow cover became 

121 automated using an algorithm that utilizes the NDSI along 

122 with a variety of threshold tests. 
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